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Abstract 

We calculate the instantaneous proper radial acceleration of test 
particles (as measured by a locally denned Lorentzian observer) in a 
Weyl spacetime, close to the horizon. As expected from the Israel the- 
orem, there appear some bifurcations with respect to the spherically 
symmetric case (Schwarzschild) which are explained in terms of the 
behaviour of the superenergy, bringing out the physical relevance of 
this quantity in the study of general relativistic systems. 
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1 Introduction 



As it is well known, since the seminal paper by Israel pQ, the only static 
and asymptotically-flat vacuum space-time possessing a regular horizon is 
the Schwarzschild solution. All the other Weyl exterior solutions [2], exhibit 
singularities in the physical components of the Riemann tensor at r = 2M. 

For not particularly intense gravitational fields and small fluctuations off 
spherical symmetry, deviations from spherical symmetry may be described 
as perturbations of the spherically symmetric exact solution 

However, such perturbative scheme will eventually fail in regions close to the 
horizon (although strictly speaking the term "horizon" refers to the spheri- 
cally symmetric case, we shall use it when considering the r = 2M surface, 
in the case of small deviations from sphericity). Indeed, as we approach the 
horizon, any finite perturbation of the Schwarzschild spacetime becomes fun- 
damentally different from the corresponding exact solution representing the 
quasi-spherical spacetime, even if the latter is characterized by parameters 
whose values are arbitrarily close to those corresponding to Schwarzschild 
metric [I]. This, of course, is just an expression of the Israel theorem (for 
observational differences between black holes and naked singularities see [5], 
[B] and references therein). 

Therefore, for strong gravitational fields, no matter how small the multipole 
moments of the source are (those higher than monopole), there exists a bi- 
furcation between the perturbed Schwarzschild metric and all the other Weyl 
metrics (in the case of gravitational perturbations). 

Examples of such a bifurcation have been brought out in the study of the 
trajectories of test particles in the 7 spacetime [7], and in the M-Q spacetime 
[8] , [9] , for orbits close to 2M [TO] , p] . 

Also, the influence of the quadrupole moment on the motion of test particles 
within the context of Erez-Rosen metric p2] has been investigated by many 
authors (see [13J and references therein). 

The purpose of this paper is to explain the bifurcation mentioned above, in 
terms of the behaviour of super-energy [H] in a neighborhood of the horizon. 
This quantity, which may be defined from the Bel [U] or the Bel-Robinson 
tensor [16] (they both coincide in vacuum), has been shown to be very useful 
when it comes to explaining a number of phenomena in the context of general 
relativity. 
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Thus, for instance, it helps to explain the occurrence of vorticity in both 
radiative [T7|, and stationary spacetimes [TS]. Also, it renders intelligible the 
behaviour of test particles moving in circles around the symmetry axis in an 
Einstein-Rosen spacetime [19]. 

In this paper we shall see how the behaviour of the instantaneous radial 
acceleration of a test particle (as measured by a locally defined Lorentzian 
observer) in a specific spacetime of the Weyl family and in regions close 
to the horizon, becomes intelligible when contrasted with the corresponding 
behaviour of superenergy. 

The Weyl metric to be considered here is the M-Q spacetime. The ratio- 
nale for this choice is that due to its relativistic multipole structure, the 
M-Q solution (more exactly, a sub-class of this solution M-Q^\ [8]) may 
be interpreted as a quadrupole correction to the Schwarzschild space-time, 
and therefore represents a good candidate among known Weyl solutions, to 
describe small deviations from spherical symmetry. 

For this metric we shall calculate the instantaneous radial acceleration of a 
test particle and the superenergy. Then the very peculiar behaviour of the 
former (close to the horizon) will be be explained in terms of the behaviour 
of the latter. 

The paper is structured as follows: in section 2 we present the Weyl family of 
metrics we shall be concerned with and briefly discuss some of its properties; 
next, in section 3, we calculate the radial four-acceleration of test particles 
in such setup. In section 4 we review and discuss the concept of superenergy. 
In section 5 we briefly describe the M-Q solution and particularize the ex- 
pressions for the proper radial four-acceleration and the superenergy for the 
case of the M-Q* 1 ) metric. Finally results are discussed in last section. 

2 The Weyl metrics 

Static axisymmetric solutions to Einstein's equations are given by the Weyl 
metric [2] 

ds 2 = e 2 *dt 2 - e" 2 VW + dz 2 ) + p 2 d<j ) % (1) 

For vacuum spacetimes, Einstein's Field Equations imply for the metric func- 
tions 

^ iPP + p- 1 ^ iP + ^ 22 = 0, (2) 



3 



r tP = p(*l-*l); 1^ = 2^*,,. (3) 



and 

Notice that (j^J) is just the Laplace equation for \l/ (in 2-dimensional Euclidean 
space); furthermore, it is precisely the integrability condition for fl3]), that is: 
given a function T satisfying ([3]) always exists. Since in the weak field 
limit ^ is related to the Newtonian gravitational potential, this result may 
be stated as saying that for any "Newtonian" potential there always exists a 
specific Weyl metric, a well known result. 

An interesting way of writing the general solution of ([21 [3]) was obtained 
by Erez- Rosen [12J and Quevedo [20J, using prolate spheroidal coordinates, 



which are defined as follows 



r _|_ r _ r _ r 



x = — ~ > 1/ 



2(7 ' 2(7 

r ± = [p 2 + (^±(7) 2 ] 1/2 (4) 

a; > 1 , -1 < y < 1, 

where a is an arbitrary constant which will be identified later with the 
Schwarzschild's mass. The prolate coordinate x represents a radial coor- 
dinate, whereas the other coordinate, y represents the cosine function of the 
polar angle. 

In these prolate spheroidal coordinates, \l/ takes the form 

oo 
n=0 

where P n (x) and Q n (y) are the Legendre functions of first and second kind re- 
spectively, and q n a set of arbitrary constants. The corresponding expression 
for the function T, may be found in |20j . 



3 The radial acceleration of test particles 

In order to find an expression for the instantaneous radial acceleration of test 
particles, it is useful to start from the geodesic equations. 

These can be derived from the Lagrangian 

2C = g a px a x^ (6) 
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where the dot denotes differentiation with respect to an affine parameter 
s, which for timelike geodesies coincides with the proper time. Then, from 
Euler-Lagrange equations it follows, 

d fdC\ dC Qj (7) 



ds \dx a J dx a 

we shall not need the full set of geodesic equations, therefore we shall display 
only the one involving radial acceleration 

2rg rr + 2r(rg rryT + g rrfi 9) - i 2 g tt ^ r - r 2 g rr ,r - 9 2 g e e,r - 4> 2 'g<M>,r = 0, (8) 

where, instead of cylindrical coordinates (p, z), we found useful to work with 
Erez- Rosen coordinates (r, 9) given by: 

z = (r — M) cos 9 

(9) 

p={r 2 - 2Mr) 1/2 sin9 

which are related to prolate coordinates, by 

r 

X= M~ l (1Q) 
y = cos 9 

Since we are concerned only with timelike geodesies, the range of our coor- 
dinates is: 

oo>t>0 r>2M n>9>0 2vr>0>O. 

Let us now consider the motion of a test particle along a radial geodesic, 
for an arbitrary value of 9. Thus putting 9 = <p — in (jSJ, and using the 
constraint (for radial geodesies) 

1 = gut 2 + g rr r 2 , (11) 

we obtain 

2rg tt g rr + r 2 (g rr g u ) >r - g u>r = 0, (12) 

It should be kept in mind that we are not interested in a full description 
of the motion of test particles (we are not going to integrate the full set of 
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geodesic equations), but only in the expression for the radial acceleration of 
a test particle at any given time. Accordingly we shall not need to take into 
consideration the constraints impossed on r, 9 and so on, from the the other 
geodesic equations, for the radial motion. 

In order to express our results in terms of physically meaningful quanti- 
ties, let us introduce (locally defined) coordinates associated with a locally 
Minkowskian observer (or alternatively, a tetrad field associated with such a 
Minkowskian observer). Thus, let 



dX = v '—Qrrdr 



(13) 



and 

It then follows that 



and 



d 2 X 

It 2 



rV-g~r 



dT = s/g^dt, 



1 - 



dX 
~dT 



dX 
dT 



(f) 2 "l 



1- (^) 2 
1 \dT ) 



dXV9rr, 

dTj ~2{-g rr f/ 2 



In the spherically symmetric case, ffTB"]) reduces to 



(14) 
(15) 



(16) 



d 2 X 
dT 2 



M 



dX\ 

It) 



2M 



-1/2 



(17) 



or, introducing the variable R = — 



d 2 X 
dT 2 



MR 3 / 2 



1 - 



dX 
~dT 



(R-2) 



-1/2 



(18) 



which is a known result. Since dX /dT is always smaller than one, the attrac- 
tive nature of gravity for any value of r (larger than 2M) is clearly exhibited 
in m 
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4 Superenergy 



As it is known, in classical field theory, energy is a quantity defined in terms 
of potentials and their first derivatives. In General Relativity however, it is 
impossible to construct a tensor expressed only through the metric and their 
first derivatives (the equivalence principle). Accordingly, a local description 
of gravitational energy in terms of true invariants (tensors of any rank) is 
not possible within the context of the theory. 

Thus, one is left with the following three alternatives: 

• Looking for a non-local definition of energy 

• Finding a definition based on pseudo-tensors 

• Resorting to a succedaneous definition, e.g.: superenergy. 

In this work we are going to explore the last alternative. As indicated in the 
Introduction, the motivations for doing so are given by the rich and profound 
physical meaning of such quantity. 

Superenergy W may be defined from either the Bel or the Bel-Robinson 
tensor [21] . Since we are working with vacuum spacetimes both definitions 
coincide, and one then has: 



W = E al3 E a p + B af3 B af3 (19) 

with 

E af3 = C^psu^u 5 (20) 

B a/3 =* C^sU^U 5 = -T) al epC ep pgU'u 6 , (21) 

where C ai ps is the Weyl tensor, r\ a ^ is the Levi-Civita tensor and u a is the 
four-velocity of observers at rest in the frame of (JTJ, i.e. 

u a = (_L 0,0,0) . (22) 



9oo 

Observe that since we are working with static spacetimes, the magnetic part 
of the Weyl tensor (B a p) vanishes identically. 

Let us next briefly introduce the metric we shall consider here, and calculate 
the corresponding expressions for the radial acceleration of a test particle 
and the superenergy. 
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5 The Monopole— Quadrupole solution, M — Q 



In [8], [9] it was shown that it is possible to find a metric of the Weyl family, 
such that the resulting solution possesses only monopole and quadrupole 
moments (in the Geroch sense [22] )■ The obtained solution (M-Q) may be 
written as follows: 



M-Q 



q" i 



(23) 



a=0 



where the zeroth order corresponds to the Schwarzschild solution. 

00 ^2n+l 



* 9 o = -£— P 2 „(cos ( 



(24) 



ra=0 



with A = Mj r and it appears that each power in q adds a quadrupole correc- 
tion to the spherically symmetric solution. Now, it should be observed that 
due to the linearity of Laplace equation, these corrections give rise to a series 
of exact solutions. In other words, the power series of q may be cut at any 
order, and the partial summation, up to that order, gives an exact solution 
representing a quadrupolar correction to the Schwarzschild solution. 

Since we are interested in slight deviations from spherical symmetry, we shall 
consider the M-Q solution, only up to the first order in q (M-Q^); with 
q > (q < 0) corresponding to an oblate (prolate) source. 

Thus, the explicit solution up to the first order, describing a quadrupolar 
correction to the monopole (Schwarzschild solution), may be interpreted as 
the gravitational field outside a quasi-spherical source, and it is given by 
(note a missprint in equation (13) in [5]) 



2 \x+ 1 



+ -g(3y 2 -l)x 



3x 2 - 1 



3y 2 -l 



In 



x — 1 

X + 1 

2x 



(25) 



3x 

{x 2 -y 2 ){3y 2 ^1) + T 
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rff-Q = I> + qT q i + q 2 T a 2 - ( I 



— qx{l-y 
8 x 2 + 1 



1 

2 

15 

32 



225 
~24 



q z In 



x- 



x 2 — y 2 



1 - ttk ( x2 + 7 v 2 - 9x V + 1 



+ 



3 x 2 — y 2 
225 



In 



x + 1 



1024 
15 



g 2 (x 2 - 1)(1 - y 2 ){x 2 + y 2 - 9x 2 y 2 - 1) In 2 
15 



x 



x + 1 



(26) 
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9(1 - y 



1 - + % 2 - 9x V + 4) 

64 



75 2 2 1 - y 2 5 2 2 1 - y 2 

— q x — n — i\ x + y 

16 x 2 — y 1 4 

^5 2/o 6 4,o42 c 2 2 i /( 2 4 4 6\ 

(2x — sc + 3s jr — orj/ + Ax y — y — y ^ 



(x 2 — y 2 ) 2 



(x 2 - y 2 f 



In [TT] it was shown that the behaviour of test particles in the M-Q*- 1 ^ space- 
time, becomes particularly strange on the symmetry axis; i.e.: 6 = 0, 7T, or 
else y — ±1 (close to the horizon). Therefore it is for that region that we are 
going to calculate the proper radial acceleration of a particle on the axis for 
the M-Q'- 1 '' spacetime (for all other regions, including the equatorial plane, 
the abnormal behaviour commented below is not observed [TT]). Using ( ITfil) 
we obtain 



d 2 X(T) _ 1 (5/4) 9j4 ( fl ) 



dT 2 



( dX(T) \ f dX(T) _ \ / (it! -2 



15gln 



V dT 
R-2 



V dT 



R 5 



R 



(R 5 - 5R 4 + 8R 3 - 4R 2 ) + 



(27) 



where 



+ q (30R 4 - 120R 3 + 130R 2 - 20R + 20) + 8R 2 - 16R 
A{R) = R(R — 2) ~ 18R2 + 8R + A+ 



1 

M 



R J 



(28) 
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Since we are interested in the value of the radial acceleration for Lorentzian 

observers instantaneous at rest, we shall plot (I27p with q = ±0.01 and 
dX(T) _ n 
dT u - 

For this metric, the expression of the superenergy (again, on the symmetry 
axis y 2 = 1) reads 

W MQ = — e (5/4)qA{R) x (768R 2 - 1152i? 3 + 576i? 4 - 96R 5 ) 
1536 

+ q (-1920/? + 3360.R 2 - 6720i? 3 + 5880i? 4 

-1200i? 5 - 540i? 6 + 180i? 7 ) 
+ q 2 (400 - 800i? + 5600i? 2 - lOOOOi? 3 + 22900i? 4 

- 32400.R 5 + 22200.R 6 - 7200.R 7 + 900.R 8 ) 

+ gin f— J (2880i? 3 - 5760i? 4 + 3600i? 5 

- 360i? 6 - 360i? 7 + 90i? 8 ) 

+ q 2 In (-2400.R 2 + 7200i? 3 - 23400.R 4 + 49200.R 5 

- 52500.R 6 + 291007? 7 - 8100i? 8 + 900i? 9 ) 

(R — 2\ 2 
— — J (3600i? 4 - 14400.R 5 + 23400i? 6 

- 19800.R 7 + 9225# 8 - 2250.R 9 + 225.R 10 ) 



1 



M 4 ( J R-2)6 J R 10 



(29) 
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Figures 
figures 



and (b) show the behaviour of W and d for q 
[c) and (d) display their behaviour for q = —0.01. 



= 0.01, whereas 
We shall next 



discuss the results obtained so far. 
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6 Discussion 



The first conclusion which emerges from figures (a) and (b) is that, close to 
the horizon, the behaviour of the test particle (in what concerns d 2 X(T)/dT 2 ) 
is extremely abnormal, as expected from Israel's theorem. 

Thus for the M-Q^ metric with q > it appears that a test particle placed 
on the axis of symmetry in the neighbourhood of the horizon does not feel 
any attraction from the source (d 2 X (T) / dT 2 ph 0). Still more shocking: as 
we move outwards (always on the symmetry axis), the magnitude of d 
increases with R, until some value of R, from which it starts to "behave" 



properly (figure (b)). 

This pathological behaviour of d 2 X(T)/dT 2 is fully consistent with that of 



superenergy in the same range of R, as indicated in figure (a) Indeed, W 
also vanishes close to the horizon, increasing as we move outwards along 
the symmetry axis, until we are far away enough from the horizon and the 
expected behaviour is recovered. 

For q < the situation is still more unusual. Indeed, on a neighbourhood 
of the horizon, on the axis of symmetry, d 2 X(T)/dT 2 > 0, implying that 
the particle experiences a repulsive force. This effect is restricted to values 
of R very close to 2. As we move away from the horizon the proper accel- 
eration becomes negative, although still displaying an abnormal behaviour 
since it increases in magnitude with R. Moving further away from the origin 
(along the symmetry axis) we recover the "normal" behaviour (d 2 X(T)/dT 2 
(negative and decreasing with R). The dependence of W with R in this 



case, displayed in figure (c) is consistent with the graphics of d 2 X(T)/dT 2 
above. Indeed, in a neighborhood of the horizon, W is singular and so is its 
derivative with respect to R, this explaining the pathological behaviour of 
d 2 X(T)/dT 2 in that range of values of R. As we move sufficiently far away 
from R = 2 we recover the expected behaviour. 

Thus we have seen that the concept of superenergy is a suitable measure of the 
strength of gravitational interaction, even in highly pathological situations. 
The fact that it is a true scalar (obtained from a true tensor) reinforces 
further its relevance in the study of self-gravitating systems. 
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